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Abstract. Magnetic field is playing an important role at all stages of star evolution from star 
formation to the endpoints. The main effects are briefly reviewed. We also show that O-type 
stars have large convective envelopes, where convective dynamo could work. There, fields in 
magnetostatic balance have intensities of the order of 100 G. 

A few OB stars with strong polar fields (Henrichs et al. 2003a) show large N-enhancements 
indicating a strong internal mixing. We suggest that the meridional circulation enhanced by an 
internal rotation law close to uniform in these magnetic stars is responsible for the observed 
mixing. Thus, it is not the magnetic field itself which makes the mixing, but the strong thermal 
instability associated to solid body rotation. 

A critical question for evolution is whether a dynamo is at work in radiative zones of rotating 
stars. The Tayler-Spruit (TS) dynamo is the best candidate. We derive some basic relations 
for dynamos in radiative layers. Evolutionary models with TS dynamo show important effects: 
internal rotation coupling and enhanced mixing, all model outputs being affected. 
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1. Introduction 

The magnetic field of a star is, like the scent of a flower, subtle and invisible, but 
it plays an essential role in evolution. Magnetic field is often not accounted for in star 
models. However, the examples below show that from star formation to the endpoints as 
compact objects, magnetic fields and rotation strongly influence the course of evolution 
and all model outputs. 

In Sect. [21 we give an overview where in evolution the fields are intervening. In Sect. 
02 we emphasize some critical observations. In Sect. SJ we focus on the general dynamo 
equations, with examples in Sect. 5 and 6 for the Tayler-Spruit (TS) dynamo. 



2. Overview on the magnetic field in star formation and evolution 

Fig. [1] shows the evolutionary track of the Sun from its formation to its endpoint with 
indications of the various effects of the magnetic field coming into play. 

— 1. Collapse and ambipolar diffusion: magnetic field may contribute to cloud 
support. For contraction to occur it is necessary that the energy density of magnetic field 
ub = B 2 /(8 it) is smaller than the density of gravitational energy uq = | R pr^g) = 

2§~rr G R l • This defines a critical mass Mb above which gravitation dominates, 

/ 5 \i $ $ 



with the magnetic flux $ = irBR 2 . In the original derivation (Mouschovias & Spitzer 
1976), a numerical factor 0.13 was obtained instead of 0.17 in this simple derivation. If 
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Figure 1. The evolutionary track of the Sun from the protostellar phase to the phase of plane- 
tary nebulae (courtesy from C. Charbonnel) with superposed indications of the various magnetic 
intervening processes 



M > Mb, large clusters or associations form. If M < Mb, no contraction occurs until 
the small ionized fraction (~ 10 -7 ) to which the field is attached has diffused (in about 
10 7 yr) out from the essentially neutral gas forming the cloud. 

2. Bipolar outflows: massive molecular outflows are often detected in region of 
star formation. A large fraction of the infalling material is not accreted by the central 
object but it is ejected in the polar directions. In massive stars, the ejection cones are 
relatively broad. Radiative heating and magnetic field are likely driving the outflows. 
Remarkably the mass outflow rates correlate with the luminosities of the central objects 
over 6 decades in luminosity, from about 1 L Q to 10 6 L Q , as shown by Churchwell (1998) 
and Henning et al. (2000). 

—3. Disk locking: from the dense molecular clouds to the present Sun, the specific 
angular momentum decreases by ~ 10 6 . Among the processes reducing the angular mo- 
mentum in stars, disk locking is a major one (Hartmann 1998). Fields of ~ 1 kG are 
sufficient for the coupling between the star and a large accretion disk. The contracting 
star is bound to its disk and it keeps the same angular velocity during contraction, thus 
losing a lot of angular momentum. The typical disk lifetime is a few 10 6 yr. 

4. Convective dynamos and magnetic braking: solar types stars have external 
convective zones which produce a dynamo. The resulting magnetic field creates a strong 
coupling between the star and the solar wind, which leads to losses of angular momentum. 
The relation expressing these losses as a function of the stellar parameters have been 
developed by Kawaler (1988). Further improvements to account for saturation effects 
and mass dependence have been brought (Krishnamurti et al. 1997). 
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Massive OB stars also have significant external convective zones which may represent 
up to 15% of the radius. Surprisingly rotation enhances these convective zones (Macdcr 
et al. 2008), which may also produce magnetic braking. 

—5. Dynamo in radiative zones: is there a dynamo working in internal radiative 
zones? This is the biggest question concerning magnetic field and stellar evolution, with 
far reaching consequences concerning mixing of the chemical elements and losses of the 
angular momentum. This question is also essential regarding the rotation periods of 
pulsars and the origin of GRBs. We devote Sect. 4 to 6 to this question. 

6. Magnetic field in AGB stars, planetary nebulae and final stages: red 
giants, AGB stars and supergiants have convective envelopes and thus dynamos. Evi- 
dences of magnetic fields up to kG in some central stars of planetary nebulae are given 
(see Jordan, this meeting), they contribute to shaping the nebulae (see Blackman, this 
meeting). White dwarfs have magnetic fields from about 10 4 up to 10 9 G, the highest 
fields likely resulting from common envelope effects in cataclysmic variables. 

3. Magnetic fields and abundances: the Henrichs et al. results 

Since OB stars also have convective envelopes, the question arises what are the possible 
fields created by the associated dynamos. If one considers a flux tube in magnetostatic bal- 
ance in the stellar atmosphere, the equilibrium field is given by the condition B 2 / (8 7r) = 
fext — -Pint- At optical depth t = 2/3, the pressure is P(t = 2/3) » (2/3) g/n. Since 
-Pint > 0, the maximum possible field B cq for magnetic equilibrium is (Safier 1999) 

B eq (r = 2/3), (^£) 1/2 - (3-D 
Table 1. The equilibrium fields. The stars are on the ZAMS, except the Sun. 

Spectral type field Spectral type field 



MO 2.8 kG GO 1.0 kG 

KO 1.5 KG F2 0.6 kG 

Sun 1.3 kG *Q9 0.2 kG 



* from the author 

Table Q] gives the corresponding estimates for stars of various types. The observed field 
intensities are often within a factor of 2 from the maximum values given in the table. 
Searches for magnetic fields in OB-type stars show no general evidence of fields above 
the level of ~100 G (Mathys 2004). This is of the order of magnitude of the possible 
fields in magnetostatic equilibrium in the convective envelopes of OB stars. The strong 
fields of 1 kG or more are not widespread (Hubrig et al. 2008). 

Noticeable exceptions were found and studied by Henrichs and colleagues (Henrichs et 
al. 2003a, Henrichs et al. 2003b). Their remarkable finding is that the few stars with high 
polar fields B p also show N and He enhancements together with C and O depletions, 
in particular for the 4 stars listed below. The abundances are given in Table [5] as the 
difference in log between the observed N abundance and the solar values. These are 
typical signatures of CNO processing, which give strong evidences of internal mixing in 
stars with a high magnetic field. These few results bring a lot of interesting questions. 

The law of isorotation of Ferraro clearly implies that an internal polar field enforces 
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Table 2. Stars with intense fields and N enrichments from (Henrichs, Neiner & Geers 2003a) 

P Ccp BUVe vsini= 27 kms -1 B p =360 G AlogN=1.2 

V2052 Oph BlIVo vsini= 63 kms -1 B p =250 G AlogN=1.3 

C Cas B2IV vsini= 17 kms" 1 B p =340 G A log N=2.6 

w Ori B2IVe v sin i = 172 kms -1 B p =530 G AlogN=1.8 

solid body rotation, see also Sect. [SI Now, the above results suggest that even in presence 
of uniform rotation, there is an efficient mixing. What is the mixing process in stars with 
a polar magnetic field? Shear turbulence is generally the main mixing process of chemical 
elements. However, it is absent here, since the stars rotate uniformly. 

The only process among those usually acting in massive stars is meridional circulation. 
Is it sufficient to produce such a mixing? In differentially rotating stars, it is usually much 
less efficient than shear mixing for the transport of the chemical elements (Meynet & 
Maeder 2000), while it is very efficient for the transport of angular momentum. However, 
we found that meridional circulation is strongly enhanced by solid body rotation, since 
uniform rotation creates a strong breakdown of radiative equilibrium. 

In evolutionary models with magnetic field and meridional circulation, there is a strong 
interplay between meridional circulation and magnetic field (Maeder & Meynet 2005): 

• Differential rotation creates the magnetic field. 

• Magnetic field tends to suppress differential rotation. 

• A rotation close to uniform strongly enhances meridional circulation. 

• Meridional circulation increases differential rotation and produces mixing. 

• Differential rotation feeds the dynamo and magnetic field (the loop is closed). 

As a result, the star reaches an equilibrium rotation law close to uniform (see models in 
the above ref.), with always a strong thermal instability amplifying meridional circulation 
and thus chemical mixing. Models show that the high enrichments in magnetic stars with 
a dynamo (Sect. [3} are essentially due to the transport by meridional circulation. Thus, 
it is not the magnetic field itself which makes the mixing, but the thermal instability 
associated to the solid rotation created by the field. 

The high surface magnetic field of these stars, which likely have a significant mass loss, 
produces a strong magnetic braking, which implies that these stars will reach a rather 
low rotation velocities during their evolution. The braking would tend to produce some 
internal differential rotation. However, the magnetic coupling is certainly strong enough 
to maintain a rotation law close to uniform, as illustrated by the models. 

4. Dynamos in radiative layers: general properties and equations 

The major question concerning magnetic fields and stellar evolution is whether a dy- 
namo operates in radiative zones of differentially rotating stars. A magnetic field has 
great consequences on the evolution of the rotation velocity by exerting an efficient 
torque able to impose a nearly uniform rotation. This influences all the model outputs 
(lifetimes, chemical abundances, tracks, chemical yields, supernova types) as well as the 
rotation in the final stages, white dwarfs, neutron stars or black holes. 

Here we first examine some general equations implied by any dynamo. The particular 
properties of the Tayler-Spruit (TS) dynamo have been studied by Spruit (2002) and we 
are using many of the equations he derived. Spruit considered the radiative zones in two 
cases, -1) when the /^-gradient dominates, and -2) when the /z-gradient is negligible. The 
more general equations of the TS dynamo have been developed by (Maeder & Meynet 
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2005). The TS dynamo is at present a debated subject. Some numerical simulations by 
Braithwaite (2006) and by Brun et al. (2007) confirm the existence of Tayler's instability. 
Braithwaite also finds the existence of a dynamo loop in agreement with Spruit's analyt- 
ical developments. However, Zahn et al. do not find the dynamo loop proposed by Spruit 
and question what may close the loop. 

4.1. Energy conservation 

If a dynamo is working in a differentially rotating radiative zone, it is governed by some 
general relations expressing the order of magnitude of its various properties. First, the 
rate of magnetic energy production Wb per unit of time and volume must be equal to 
the rate W v of dissipation of rotational energy by the magnetic viscosity v. We assume 
here that the whole energy dissipated is converted into magnetic energy. The differential 
motions are those of the shellular rotation with an angular velocity f2(r), so that the 
velocity difference at radius r is dv — r df2. The amount of energy corresponding to a 
velocity difference dv during a time dt for an element of matter dm in a volume dV is 

W v = ±dm(dv) 2 ^j t = ± Q v(^j =\qvVl\ 2 with g = r|Vn|/n, (4.1) 

because the viscous time dt over dr is given by dt = (dr) /v. The magnetic energy density 
is ub = B 2 /(8tt), it is produced within the characteristic growth time of the magnetic 
field cr^ 1 , thus the rate Wb of magnetic energy creation by units of volume and time is 

B 2 1 

W B = —o-b = -u\r 2 o~BQ , (4.2) 

where we have used the expression of the Alfven frequency uja — ■ Now, let us 

assume W v — Wb, i.e. that the excess of energy in the differential rotation (compared 
to an average constant rotation) is converted to magnetic energy by unit of time. This 
gives the following expression for the viscosity coefficient of magnetic coupling 

This is the coefficient which intervenes in the expression for the transport of angular mo- 
mentum, in the Lagrangian form as given by Eq. (|5.8|) below. Let us note that compared 
to the energy available for the solar dynamo driven by convection, the amount of energy 
available from differential rotation is very limited. 

4.2. The a and Lo-effiects: vertical instability and stretching of the field lines 

A dynamo needs both the a-effect and w-effect. The a-effect consists in the generation 
of a poloidal field component from the horizontal component. In the Sun, the a-effect 
is created by the convective motions and by the twisting of the magnetic loop by the 
Coriolis force. However, other instabilities with a vertical component may produce the 
necessary a-effect. The w-effect consists mainly of the stretching of a small radial field 
component in the East-West directions. The winding-up of the field lines generates a 
stronger horizontal field component, converting some kinetic energy into magnetic energy. 

If due to an instability in radiative layers, some vertical displacements (necessary for the 
a-effect) with an amplitude l r /2 occur around an average stable position, the restoring 
buoyancy force produces vertical oscillations with a frequency equal to the Brunt-Vaisala 
frequency N. The restoring oscillations will have an average density of kinetic energy 
m n = /jv P^ 2 N 2 , where /at ~ 1. In order to produce a vertical displacement, the magnetic 
field must overcome the buoyancy effect. In terms of energy densities, this is ub > 
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un, where ub has been given in the previous section. Otherwise the restoring force of 
gravity would counteract the magnetic instability at the dynamical timescale. From this 
condition, one obtains I 2 < r 2 jfe. If, /n = |, we have the condition (Spruit 2002) 

I < lr = r ^ , (4.4) 

where r is the radius at the considered level in the star. 

The stretching of the field lines for the w-effect is governed by the induction equation 

— =V x (vx B)+nV 2 B . (4.5) 

An unstable vertical displacement of size £ from the azimuthal field of lengthscale r and 
intensity B v also feeds a radial field component B r . The relative sizes of these two field 
components arc defined by the induction equation (14. 5p , which gives the following scaling 
over the time St characteristic of the unstable displacement, 

B r rs SB w — — B v St . (4.6) 
r St 

For the maximum displacement l r given by Eq. 14.41 this gives 

which provides (Spruit 2002) an estimate of the ratio of the radial to azimuthal fields. 

4.3. The magnetic and thermal diffusivities 

The magnetic diffusivity n tends to damp the instability, while the thermal diffusivity 
K produces heat losses from the unstable fluid elements and thus reduces the buoyancy 
forces opposed to the magnetic instability. Both effects have to be accounted for. 

If the radial scale of the vertical instability is small, the perturbation is quickly damped 
by the magnetic diffusivity r\ (in cm 2 s _1 ). The radial amplitude must satisfy, 

lr > — , (4-8) 

where, as seen above, <tb is the characteristic frequency for the growth of the instability. 
The combination of the two limits (|4.7p and (|4.8[) gives for the case of marginal stability, 

For given rj and ctb, this provides the minimum value of wa, and thus of the magnetic 
field B, for the instability to occur. The instability is confined within a domain, limited 
on the large side by the stable stratification (|4.4[) and on the small scales by magnetic 
diffusion (|4.8[) . For the case of marginal stability, which is likely reached in evolution, 
this equation relates the magnetic diffusivity 77 and the Alfven frequency uix- 

The Brunt -Vaisala frequency N of a fluid element displaced in a medium with account 
of both the magnetic and thermal diffusivities n and K is (Maeder & Meynet 2004) , 

N 2 = ^ro N 't.*<± + N 1> ( 4 - 10 ) 

K 

with iVf, ad = j^- (Vad — V) , and N 2 = ^V^, (4.11) 

Hp Hp 
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The ratio rj/K of the magnetic to thermal diffusivities determines the heat losses. The 
factor of 2 is determined by the geometry of the instability, a factor of 2 applies to a thin 
slab, for a spherical element a factor of 6 is appropriate (Maeder & Meynet 2005). 

4.4. The magnetic coupling and the timescale Ob 

The momentum of force S by volume unity due to the magnetic field is obtained by 
writing the momentum of the Lorentz force F^. The current density j is given by the 
Maxwell equation ^j = VxB. Thus, one has 

S = f x F L = -fx (jx B) = — fx f (V x B) x B] , (4.12) 

C 47T V / 

in modulus S « ±- B r B v = J- ( l A B% = p r 2 (^] . (4.13) 



An 1 v ~ 4tt \rj H \N , 

The units of S are g s -2 cm -1 , the same as for B 2 in the Gauss system. The kinematic 
viscosity v (in cm 2 s _1 ) for the vertical transport of angular momentum is 

77 1 dr 1 dr 1 d In r , 

u = -L = -F— = -F = -F , (4.14) 

q g dv g rdfl g il d In fl 

where F is a force by surface unity, which also corresponds to a momentum of force by 
volume unity in g s -2 cm . F is applied horizontally to a slab of velocity v in a direction 
perpendicular to r. Considering only positive quantities, with q = \dhiil/dhxr\, one has 

S *' (4.1 5 ) 



pqtt Nqfl 

Now, we can compare this expression for v to Eq. (|4.3[) and get 

o-b = N ■ (4.16) 

This important expression relates the growth rate of the magnetic field to its amplitude 
(through wa). It can also be obtained by expressing the amplification time r a of the field 
line B r to the level of B v by the winding-up of the field line B v w B r r (— — ] - 



dr . 

N 



waO q 

and equaling this to Cg 1 we also get Eq. (|4.16p . 

4.5. The basic equations 
Introducing the expression (|4.16|) of ctb in Eq. (|4.9p . we get 

2r> /^A^ 3 



(4.17) 



v = r 2 nq^) . (4.18) 



2 



Also, with the expression (|4.10[) . we can write for a 

2 uln 2 q 2 , 

a* = jr . (4.19) 

k w 2 4- A^ 2 

5+2 iv T, ad ' JV M 

These equations are quite general. If the growth rate a of the instability is known, the 
two equations (|4.18p and (|4.19[) form a system of 2 equations with 2 unknowns rj and wa- 
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5. The case of Tayler Spruit dynamo 

In a non-rotating star the growth rate of the Tayler instabilty is the Alfven frequency 
wa- In a rotating star, the instability is also present, however the characteristic growth 
rate <tb of the instability is, if wa *C f2, 

= f , (5.D 

because the growth rate of the instability is reduced by the Coriolis force (Spruit 2002). 
If so, Eqs. (|4~T8f and (|4~19|) become 



/WA\ 2 _ » 2 g 2 



V O / AT 2 r '/ K l N2 ' 

JV T,ad V /K + 2 + J > 



g 2 V n J 



6 



(5.2) 



(5.3) 



This forms a system of 2 equations for the 2 unknown quantities 77 and cja- With a new 
variable x = (cja/^) 2 , we get (Maeder & Meynet 2005) a system of degree 4, 



r 



„ > ^ K + ^) ^ " + 2^ x - 2ftV - . (5.4) 

The solution x provides the value of the Alfven frequency uia and thus of the B field. By 
(|5.3|) one gets the value of r\ and by (|4.15|) the value of v. The above equation applies 
to the general case where both JV M and Nt are different from zero and where thermal 
losses may reduce the restoring buoyancy force. The solutions of this equation have been 
discussed (Maeder & Meynet 2005). In particular, if Nx = 0, one has 

*=( 9 ^) 2 and v=r2nqi (£) 6 ' (5 - 5) 

which shows that the mixing of chemical elements decreases strongly for larger /j gradients 
and grows fast for larger q values. 

The ratio wa/O given by the solution of (|5.4|) has to be equal or larger than the 
minimum value defined by (|4.9[) . This leads to a condition on the minimum differential 
rotation for the dynamo to work (Spruit 2002), 



When iV 2 is larger, as for example when there is a significant \i gradient, the differential 
rotation necessary for the dynamo to operate must also be larger. If the above condition 
is not fulfilled, there is no stationary solution and the dynamo does not operate. In 
practice, this often occurs in the outer stellar envelope. 

5.1. Equations of transport of chemical elements and angular momentum 

The equation for the transport of chemical species with mass fractions X% is at a La- 
grangian mass coordinate M r , 
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WITH MAGNETIC FIELDS 

X c = 0.70 
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Figure 2. Left: evolution of the angular velocity Q as a function of the distance to the center 
in a 20 Mq star with Vi n i = 300 km s" 1 . X c is the hydrogen mass fraction at the center. The 
dotted line shows the profile when the He-core contracts at the end of the H-burning phase 
(Meynet & Maeder 2000). Right: rotation profiles at various stages of evolution (labeled by the 
central H content X c ) of a 15 Mq model with X = 0.705, Z = 0.02, an initial velocity of 300 
km s _1 and magnetic field from the TS dynamo (Maeder & Meynet 2005). 

where -D e ff is the coefficient for the transport by meridional circulation and the possible 
horizontal turbulence. The equation for the transport of angular momentum is 

where U2 (r) is the amplitude of the radial component of the velocity of meridional circu- 
lation and v the value given by (|4.3[) . This equation is currently applied in stellar models 
for calculating the evolution of f2. With account of the detailed expression of ^(r), which 
contains terms up to the third spatial derivative of fi(r, t), the above equation is of the 
fourth order and its numerical solution requires great care. 

6. Numerical models 

Numerical models accounting for meridional circulation and magnetic field generated 
by the TS dynamo have been computed (Maeder & Meynet 2005). The resulting fields are 
a few 10 4 G through most of the envelope, with the exception of the outer layers where 
differential rotation is too small to sustain the TS dynamo. The diffusion coefficient for 
the transport of angular momentum is large. In the Sun, it is of the order of 10 2 to 10 6 
cm 2 s -1 , sufficient to impose solid body rotation at the age of the Sun (Eggenberger et 
al. 2005). This coefficient is much larger in more massive stars, in the range of 10 10 to 
10 12 cm 2 s _1 in a 15 M Q star. There, it imposes nearly solid body rotation during most 
of the MS phase, while without the field there is a high differential rotation (Fig. [5]). 

The nearly solid body rotation of star with magnetic field drives meridional circulation 
currents which are faster than the currents in differentially rotating stars. This leads to 
large surface enrichments in N and He together with C,0 depletions in massive stars 
(Fig. [3|). Thus, the enhanced mixing results from the thermal instability enhanced by 
uniform rotation. The stellar lifetimes are enlarged by the mixing and the other model 
outputs are also modified (Maeder & Meynet 2005). Therefore, magnetic field is also a 
basic ingredient of stellar evolution. 
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Figure 3. Time evolution of the surface helium content Y s in mass fraction, of the N/O, N/H and 
N/C in mass fraction for various models: The dotted line applies to the model without rotation, 
the short-broken line to the model with rotation (vini= 300 km s _1 ) but without magnetic 
fields, the continuous line to the model with rotation (vini= 300 km s _1 ) and magnetic fields 
from the TS dynamo (Maeder & Meynet 2005). 
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